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ABSTRACT 


Recent improvements in rare earth magnets have made it 
fossible to construct strong, lightweight, high horsepower 
DC motors. This has cccasioned a reassessment of electrome- 
chanical actuators as alternatives to comparable pneuratic 
and hydraulic systems for use in flight control actuators 
for tactical missiles. This thesis cevelops a lew-crder 
mathematical model for the simulation and analysis of brush- 
less DC motor performance. The model is implemented in CSHP 
language. It is used to predict such motor performance 
curves as speed, current and power versus torgue. 
Electronic commutation based on Hall effect senscr_ rosi- 
tional feedtack is simulated. Steady state motor kehavior 
is studied under both constant and variable air gap flux 
conditions. The variable flux takes two different forms. 
In the first case, the flux is varied as a simple sinusoid. 


In the second case, the flux is varied as the sum of a sinu- 
soid and one of its harmonics. 
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I. INTRODUCTION 


Direct current electric motors are jaining wider use in 
space aprlications where long life, hign reliability, high 
torgue and liyht weight are critical factors. In recent 
years, advances in magnet technology ani materials, esfye- 
cially the use of rare earth magnets, have made fossitle 
significant imprcvements in the torgue to inertia ratios of 
permanent magnet (PM) motors. Samarium cobalt is cne such 
rare €arth magnetic material. The great potential of the 
Tare e€arth magnets derive from their inherent high flux 
density and high coercivity properties. Higher flux densi- 
ties mean greater developed motor torgue while nigh coer- 
civity means greater innate resistance to demagnetization 
which permits thinner magnet sections. The two factors in 
combination result in increased mechanical power and energy 
product with decreased physical size and weicht. 

Brushless dc motors belong to this class of permanent 
magnet motors and thus enjoy their improved torgue to size 
Characteristics. In fact, brushless dc motors are sigrifi- 
cantly smaller than either their ac counterparts or conven- 
tional brush-type dc motors on an egual horsepower basis. 
Compared to its conventional counterpart, a brushless ic 
motor*s structure is “inside out." In the brushless config- 
uration, the permanent magnets are attached to the rotecr and 
the ccnducting coils are placed in the stator. As its name 
suggests, brushes have been eliminated from this tyre of 
motor. Instead, ccmmutation of current in the stationary 
armature is accomplished electronically by switching cr the 
appropriate windings via transistors specially designed to 
provide the high currents needed for motor pover. In crder 
to seguence the current to the coils, the commutation 
circuitry is provided rotor position feedback from sensing 
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devices, most commonly from Hall effect sensors. These are 
simple semiconductor devices that develop a jpeclarizei BE 
voltage depending ufon the magnetic field passing through 
them. Small, highly sensitive and reliable, Hall effect 
senscrs reguire little power to operate and thus are widely 
used. cai 
Erushless dc motors possess ali the advantages tradi- 
tionally associated with conventional dc motors such as : 
linear torque-speed characteristics, excellent response 
times and high efficiency. In addition, there are a nuacer : 
af of important advantases that brushless dc motors pcssess 
over standard brush-type dc motors. An obvious advantage is 
that electronic commutation means the elimination of comnmu- 
tators, the wear and tear of brushes, and the build up of : ; 
carbon dust which heretofore significantly limited the life aa. 
of the standard dc motor. Since electronic current 
switching essentially involves no component wear, there 
7 results improved reliability and enhanced lifespan. As weil 
~~ as being subject to wear ani tear, the action of brushes 
. sliding over commutators also generates radio frequency 
interference (RFI) which is a severe liability in many 
* applications. In fact, it is often in hostile and exrlosive 
i environments that tke advantages of dc motors are most 
needed. Another benefit is in heat transfer charactezis- 
tics. By placing windings in the stator the thermal pati to 
the environment is made shorter and more direct. with the 
removal of heat thus enhanced, mechanical and electrical 
malfunctions due to keat are reduced. A further important 
advantage of the brushless dc motor is in the area of motor 
= speed variability. Whereas ac motors usually operate at one 
speed fixed by the powerline frequency, azy dc motor's speed 
can ke ckanged to meet varying power requirenents by singly 
adjusting the dc input voltage. The electronic commutation 


scheme cf the brushless de motor lends itself nicely to just ae 
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this type of variabie speed control. It is a relativel; 
Simple matter to attach to the commutation circuit tEcari a 
snall, inexpensive microprocessor whose logic can be 
programmed to run tke motor at any of several speeds or in 
either direction. A single motor can thus be proyranned, in 
place if desired, so that its speed, direction and tergue 
are matcted to a particular application. 

There are at least aS many applications for frushiess 
dc motors as there are for brush-type motor systems as weil 
as a hest of new apflications. Of particular relevance t? 
this thesis is the use of brushless dc motors in <Zlaight 
contrel actuator sytems of Navy tactical missiles, in 
particular, the cruise nissile. There are three gereral 
types of actuators that are in use in missile systems: fneu- 
matic, hydraulic and electromechanical. The maneuvering 
requirements of the mere advanced tactical missiles call for 
high contrel torques which heretofore could only be produce? 
by high pressure pneumatic or hydraulic actuators. Bcth of 
these actuator systems are more susceptible to mecnanical 
malfunction than their electrical counterparts. With the 
improved torgue to inertia ratios resulting from better rare 
earth magnet materials and the higher mechanical reliability 
due to removal of brushes, as well as other advantajes 
previously enumerated, the brushless de motor is ina fosi- 
tion to compete with most pneumatic and hydraulic systems. 
This thesis is but one part of a detailed study of the state 
of the art in electremechanical actuators as applied te the 
cruise missile system (Ref. 1]. 


A. FEROGEAMNING LANGUAGE 


The objective of this thesis is to create a simplifiei 


low-crder mathematical model that will 


the response of a brushless de motor. 


accurately simulate 
The IBM Continuous 


System Mcdeling Program (CSMP) was chosen for use since it 


is a user-oriented ccmputer language 


Specifically created 


for modeling dynamic fhysical systems and is one of the mest 


widely used in this area [Ref. 2]. 


It has the flexibility 


or determining the response of a system that is w2odeled 


either in a block diagram format or 
differential equations. 


as a set of ordinary 


While CSMP is an offshoot of IBM's 


Digital Simulation Language, it embodies Fortran as its 


source language and retains uch of Fortran's capability and 


flexibility. CSMP is an applications-oriented preyran 


intended for use with the IBM mainframe systems. the 


utility and ease of using this program 


simplified frogram ccntrol statements 
describe the mathematical equations or 
the system and from the flexibility of 
which contains preprogrammed function 


languajye stems fron 
that almost exactly 
physical variables of 
its prograg structure 
blocks that obviate 


the need for complicated Subroutine programminy. In 


essence, CSMP was chosen because it 


permits the user to 


concentrate on the details of the physical system rcather 


than on the time-consuming complexities of programming. 


EB. SYSTEM BLOCK DIAGKAS 


The block diagram format was chosen as the basis for 


the CSMP program since it is a convenient way of analyzing 


the input-output relationships cf a physical system and the 
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flow of signals through it without having to refer to the 
system equations. The basic idea of a block diagram stems 
from the application of Laplace transforms to a system's 
integro-differential cr differential eguations. In effect, 
a differential equaticn is rendered into a simpler aljetraic 
expression which, in turn, becomes a fundamental equaticn of 
a system's functional blocks. Each block indicates the 
relationship between its input driving signa and output 
response; ie., the transfer function of the input andi 
output. The transfer function is defined as the ratio of 
the Laplace transform of the output to the transform of the 
input with all initial conditions assumed to be zerc. A 
block's transfer function does not include any information 
about the internal structure of that part of a systen, 
merely the transformation of a signal between two pcints in 
a systen. Furthermcre, it should be borne in mind that all 
physical systems have certain transfer characteristics that 
are actually nonlinear to some extent. A dic motor systen 
can be fairly accurately described using the transfer func- 
tion approach to modeling and, when nonlinear transfer char- 
acteristics exist, CSMP modeling can usually be obtained 
from preprogrammed ncnlinear function routines. 


C. ODEVEICPMENT OF ACTOR SYSTEM EQUATIONS 


In crder to construct the block diagram of a systen, 
one affroach is to identify the various mathematical reia- 
tionships between the components of the system and to write 
balance eguations that will identify each individual block 
comprising the system. The system equations for a brushless 
dc motor are presented in chapter three where all the 
details of the logic control and switching of transistors 
for cemmutation of tke motor are developed. ‘hat follows is 
a Simplified input-output characteristics model of the type 
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that has 
contrel systen 
input-output model 
modeling 

additional 


mechanical engineering study cf load torgue reyuirenents tor 


the mctor [Ref. 1]. 

To write 
fasic de 
Kirchoff's voltage 
circuit 


cH ---- 


R 


Figure 2.1 


been a Standard for 
studies [Ref. 3]. 


frocess for 


usefulness of 


the electrical 
Motor, the 
law is applied 
(Figure 2.1). 


used in 
of this 


step in the 


brush-type dc motors 
The development 


the first learning 


the brushless dc motor and has) the 


being a tool for a concurrent 


balance equation describin,; a 


classical loo: method based on 


to the eyuivalent motor 


mc a a a ee 
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Equivalent DC Motor Circuit 
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The Easic dc notor has the coil windings on the armature and 
the magnetic field scurce on the stator. In the case of the 
brushless dc motor, the rotor consists of permanent nagnets 
that produce the magnetic field and the stator contains the 
coil windings. The Kirchoff Law for the stator circuit is 


es(t) = Ldi/dt + Ri(t) + eb(t) (egr 2.1) 


where R and L are the stator resistance and inductance 
respectively, and eb(t) is the back enf. When a conductor 
moves in a magnetic field, or there is any relative 2ction 
ketween the magnetic field source and the conductor, a 
voltage is generated across the terminals of the conductor. 
In the case of the dc motor, the voltage is proportional to 
the shaft velocity and tends to oppose the current flow. 
The relationship between the back emf and the shaft velccity 
is 


es (t) (Km*) *wm(t) = Kb*wm (t) (eqn 2.2) 


es (t) 


Kb*dOm (t) /dt (egn 2.3) 


Substituting and rearranging equations 2.1 and 2.2 into 
their differential forms gives the following two equations. 


disdt=(1/L) *es (t) - (F/L) *i (t) -(1/L) *Kb*um (t) (eqn 2.4) 
dem (t)/dt = wm(t) (eqn 2.5) 
The basic balance equation that describes the motor as 
a mechanical system derives from application of Newton's 
laws of motion to the system components. The dynamic eyua- 


tion for a motor coupled to its load is 


tm(t) = Jedwm(t)/dt + Bewm(t) + t1(t) (eqn 2.6) 
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where tl(t) is the load torque. 


recomes 


Rearranging the equation it 


dwm (t) /dt= (1/3) *tm (t)- (1/9) *tl (t) - (B/J) *we(t)  (egn 2.7) 


Bis the total viscous friction 
is comprised of the sum of the 
motor, Bm, and the load, Bl, 
Likewise, J is the total systen 


of the motor inertia, Jm, and 


reflected through the coupling 


coefficient of the mcter and 
viscous friction due tec the 
as seen by the motor shaft. 
inertia comprisea of the sua 
of the load inertia, Jl, 


device to the motor shaft. 


Foth Bl and Jl are functions of the motor's speed relucing 


mechanisg where N is the speed 


reduction ratio (N>1). The 


following equations summarize the foregoing. 


B = Bm + BL 


J = Jn + Jl 
BL = Blp/N 
Jl = J1p/Nn 


(eqn 2.9) 


(egn 2.10) 


(egn 2.11) 


Jlp and Blp are the inertia and viscous friction seen at the 


load side of the systen. Since a dc motor is essentially a 


torque transducer that converts electricai energy ito 


mechanical energy, the following equation is necessary to 


show the relationshif amoung the developed torgue, tz(t), 


the air gap flux, $, and the stator current, i(t). 


tm(t) = (Km*) *i(t) 


(egn 2.12) 


Because the magnetic field in the motor is assumed uniforn 


and constant, this can be simplified to 
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tm(t) = Kt*i (t) (egn 2. 13) 


where Kt is the torgue constant of the motor. 

Eguations 2.1 through 2.13 represent the cause anh 
effect eguations of the systen. The application of es (t) 
Froduces a current flow which causes the torque and the back 
emf te be generated. The torgue produced then causes the 
angular displacement Om(t). Given the differential ejua- 
tions of the system, the Laplace transform can te aprlied 
and the block diagram of the system jenerated. AS an 
exaaple, the block fer equation 2.1 is generated as follows 


es (t) - eb(t) = L¥di/sdt + R*i(t) (egn 2. 14) 


Aies(t) - eb(t)} =L(L*disvdt + R*i(t)} (egr 2.15) 


En(s) = Es(s) - Eb(s) = (Ls + R)*I(s) (egn 2. 16) 
I(s)/En(s) = Is/{Es - Eb) = 1/(Ls + R) (eqn 2.17) 


Performing the same operation to the remaining differential 
equations, the block diagram in Fiyure 2.2 resuits. 

It is important to note that while a dc motor is tasi- 
cally an open-loop system, the back emf of this tyre of dc 
motor acts as a natural feedback loop that tends to improve 
the stability of the motor. 

Cnce the basic frlock diagram is developed, writing a 
CSMP frogram is simple and fairly straightforward. Care 
must be taken, however, to put the block transfer functions 
into the proper format for CSMP*s functional blocks. A copy 
of the CSMP program fora tkasic dc motor is provided in 
Appendix A. 


f | 
. ANGULAR 
DISPLACEMENT 
a 
| 
| 
Figure 2.2 Block Diagram of DC Motor = 
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Pe A. OVERVIEW 


fi Cne of the purgeses of this model is to produce a set 
of performance curves that accurately simulate a given 
: Motor's behavior when operated under varying load cendi- 
tions. As such the model can be a useful tool for studying 
(. the changing performance due to the configuration changes of 
a motor under development or for observing the behavior of a 
: motor under varying conditions of application or environ- 
ment. Thus, having developed an input-output computer model 
for the DC motor, the next step was to assign values tc the - 
model parameters and run the progran. Program inruts, 
constants and parameters were assigned values that are 
typical of brushless DC motors commercialiy available. 
Table I rrovides the parameter values for a given motor that oe 
will te used as a generic motor for the following analysis 
and discussion. 


TABLE I 
Typical Ccmmercial Motor Parameters 


Stator Resistance, : ohns 


2.74 
Stator Inductance 0.0016 henries 
Torgue constant, kee 15.9 oz-insame 
Back EMF constant, Kb Q.112 volt/rad/s 
Rotor Inertia, Jn 0.001 oz-in/s?2 
Viscous Friction 
Coefficient, Em unknown ozZ-inysrad/s 


ies as gy ngs CE tae ge 8 ca ee eggs Se ones ete tae te 
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B. SPEED VS TORQUE CURVE 


Specifically, three curves were produced: motor speed 
vS. genéerate? torque, current vs. torgue and output fower 
vs. torque. Since permanent magnet DC motors rave Linea 
characteristics, a straight-line curve resultei1 under all 
load conditions when the parameters of Table I were inserte 
into the model and a voltage of 30 VDC was anvplied. In 
order to fully understand the individual contributions of 
the varicus motor parameters on tne speed-torcue curve, it 
was necessary to determine the mathematical relaticnshifs 
btetween the parameters. Feferring to equations 2.1 and 2.5, 
the electrical and motor eguations when considered urder 
steady state conditicns becone 


es (t) P*i(t) + Kb4wm(t) (eqn 3.1) 


tm(t) = Xt*i (t) (egn 3.2) 
Solving for i(t} and combining the equations resuits in 
es(t) = (R/Kt) *tm(t) + Kb*wm (t) (egn 3.3) 


When there is no torgue applied to the motor, the current is 
small encugh to neglect and the no-load velocity becomes 


wnl = es/Kb (egn 3.4) 


Given that the applied voltage is a constant, it is apparent 
that knowledge of tke value of the back emf constant, Kb, 
will previde the various no-load speeds for any value of 
applied voltage. On the other hand, in the absence of any 
known motor specifications, the no-load speed of a given 
motor can te empirically measured and then the back enf 
constant, Kb, can be mathematically derived. 
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Looking next at tke slope of the speedi-torgue curve, it 


can ke seen that whereas one endpoint of the curve ternmi- 
nates at no-load speed, the other ends at stall which is 
that value of torgue which is large enough to literally stop 
the motor. From equations 3.1 and 3.2, with wm set to zero 
the generated torgue, tm(t), at stall is ts given by 


ts = (Kt/R)*es {eqn 3.5) 


Combining equations 3.4 and 3.5, the equation of a straight 
line for the speed-torque curve results in 


wo(t) = wnl - (Rm*tm(t)) (egn 3.6) 
where 
Rm = R/(KbE*Kt) (eqn 3.7) 


and km is the slope factor for the curve, also called the 
speed regulation constant. Oktviously, changing any or all 
of the parameters will alter the slope factor, Rm. However, 
the Lack erf constant, Kb, and the torgue constant, Kt, are 
strictly related by the common factor, aizc gap flux, ¢. In 
fact, when both constants are put in the MKS system of 
units, Kb in volt/rad/s and Kt in Nm/amp, then they are 
egual as equation 3.8 shows. 


Kb(voltyrad/s) = Kt (Nm/amp) (eqn 3.8) 


Thus, knewledge of either constant provides knowledge of the 
ether. Referring to equation 3.7, it can be seen that the 
final unknowns are tke resistance, R, and the slope factor, 
Rm. If the resistance, R, is known from manufacturer srpeci- 
fication sheets, for example, then the slope, Rm, can be 
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easily sclved. On the other hand, if siven an empirically 


Measured speedi-torgue curve, it is now obvious that the 
curve not only provides Kb which in turn provides Kt, but it 
also can provide a derived value for motor resistance, R. 
Cnce acquired, the three basic ,araneters, Kb, Xt andi &, 
were plugged into the CSMP model and the speed-torgque curve 
for a 30 VDC input was generated (see Figure 3.1). 

This same procedure can be used to produce a family of 
speed-torqgue curves if it is desired to study a motors 
performance under different input voltages. Figure 3.2 
shows the family of curves produced when the values given in 
Table I were entered into the model. 


C. CURRENT VS TORQUE CURVE 


With the speed-tcrgue curve understood, the next task 
was to repeat the same general procedure for the frotor 
current vs. torque. As a starting point, the equation for 
the steady-state current is given by 


iss(t) = (es(t) - Kb*wm(t))/K (egn 3.9) 


Since KE and R were fixed by the speed-torgue curve, gener- 
ating a current-torgue curve was straightforward procedure 


for any applied voltage (see Figure 3.3). 


D. CUTPUT POWER VS ICRQUE CURVE 


The final performance curve to he studied was the 
output pcewer vs. torgue curve. The basic equation fer fower 
is given by 


P(t) tm(t) *wm(t) (egn 3.10) 


P(t) 


Kt*i(t) *wm (t) (egn 3.11) 
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@ Figure 3.1 Motor Speed vs Torque Curve 
* where pcwer iS measurel in watts. Since this eyuation 
oe involves no new constants, jprojiuciny the power curve «sas 


only a watter of addiny a line to the existing aodel which 
also included a conversion factor for changing the tervus 


units frem oz-in/amp to Nm/amp in order for the power to he 
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Figure 3.2 Family of Speed-Torgue Curves 


dimensionally correct in watts. Again cunning the peejran 
with the parameters given in Table I, the model }roduced the 


Fower curve given in Figure 3.4 
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Figure 3.3 Motor Current vs Torque Curve 

As expected both the speed-toryue and current-torjue 
curves were linear. The power-torjue curve is obviously not 
linear, especially at high lcais where the power Lejins to 
roll cff fairly shargly. From the molelling peint of view 
this is due solely tc the fact that the motor sicws decwna 
faster than the load torque increases {see eyuation 3.9). ae 
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Figure 3.4 Output Power vs Torque Curve 


Another factor that further accentuates this power rollcff£ 
in an actual motor is the effect due to armature reaction. 
Whenever a current flows through the motor armature ina 
permanent mayneét dc mctor, the armature Lecones an electro- 
magnet which produces a flux that tends to oppose the fiux 
Froduced by the fermanent magnets. This has the effect of 


e! 


partially jlemagnetizing the permanent maguets. Tnis 
demagnetization is a reversitle effect, meaning that as 
current returns to zero the permanent magnets returr t9 
their full strength. However, as the current increases and 
armature reaction sets in, it has the effect of decreasing 
the torgue constant, Kt (as well as Kb). AS equation 3.35 
indicates, this also contributes to the output power rolicff 
at high loads. Though current zermanent tagret CC actors 
made cf rare earth magnets have high coercivity taat resists 
armature reaction, at high levels of rated current a small 
amount of armature reaction still occurs. For the purroses 
of this thesis, the model is considered to accurately sizu- 
late motor behavior fcr loads up to the peak power load tor 
the given applied vcltage. From that load upward, it is 
assumed that armature reaction is likely to occur and result 


in nonlinear behavior which is not included in this rodel. 


Fe MCTOR REVERSAL 


The final step in the modeling frocedure was to reverse 
the motor's direction and ensure that wirror images or ths 
three curves resulted. To do this, it was first necessary 
to replace the positive input voltage with a negative ccunt- 
erpart. Next, in order to maintain mocéi consistercy where 
the lcad tcrque, tl(t), is treated as a torzgue that opposes 
the meticn of the mctor, it was necessary to modify the 
model so that the load torque was a pcesitive value. Tale 
consisted of creating a CSMP procedure block that identifie3 
the input voltage as either positive or negative ard then 
treated the load torgue accordingly. Once program tugs were 
removed, the nodel accurately simulated notor reversal wita 
proper speed, current and power values for the appropriate 
load conditions. 


IV. ELECTRONIC COMMUTATION fe 


A. SWITCHING ACTION 


The next major step in this thesis was to convert the 
: basic model of a standard brush-type dc motor into a three- 
phase, fcur-pole brushless dc motor systen. Figure 4.1 
shows this systen with the three-phase stator windings 
configured in a standard 'star' connection with each winding 
awd oriented 120 electrical degrees from tne others. 


Figure 4.1 Controller Configuration 


The six transistors are connected to the ends of each stator 
leg and through logic-controlled switching action previde 
three-phase full-wave motor ccntrol. Figure 4.2 indicates 

the switching logic sequence for counterclockwise rotation, 
'e To produce clockwise rotation, the sejuence is reversed. 
In either case, a pair of transistors will be switched at 
the start of each 30 degree (mechanical) interval which 


causes current to Sisultaneously build up in one leg, flow 
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Figure 4.2 Switching Logic of a 3-Phase ®rushless CC Motor 


steadily in a second, and decay to zero in a third. Given 
proper sequencing, the net deveioped torjie ideally ceaches 
a steady state valte that causes motor rotation at a 
constant speed in the desired direction. 

For windings configured in a star arranjement, it can 
be seen that conduction is ccntinuous in one leg while 
commutation occurs in the other two lejs. Tor exauyir, as 
Q1 and @5 are energized between 30 ani 60 idejrees (nechan- 
ical), current flows down leg A ani into leg Be. At the sage 
time, the current through leg C jue to the energy stcrag» 
behavior of an inductor decays to zero throwjh dicie, 3? 
(see Figure 4.1). In the next seguence, 60 to 90 deygrercs, 
C6 is energized and ¢S switched off. Current has reach 1 a 
steady state flow threugh leg A, but now it flows dewn tej 7 
and the current in leg B decays to zero through D2. In the 
next sequence, 60 tc 120 degrees, Q2 is Switched on aid ¢1 
Switched off. Current has reached steady state flow through 
leg C but now current builds up and flows fron leg B wiile 


the current in leg A decays through lej C. This same action 
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Tepeats in subseguent commutation states. Lookiny at roint 


Cin Figure 4.1, it can be seen that, Since current is 
always continuous in one leg whether flowing into cr f£rozr 
node C, Kirchoff's Current Law requires tzat the net current 
flowing in the other two legs must eygual the flow in the 
continuous leg. This permits the following inperta 
Simplification which was needed later in the tofeling 
Frocess. Namely, the three-legged stator can be approx 
Mated by just two windings serially configured. Ore le 
wili always accurately reflect steady state flow conjiticns 
while the other will be treated as if it were in ste 
state because the build up in one leg is balanced ry deca: 
in the cther. It is recognized that this nodel is a first 
approximation and that future modeles will regnire adii- 
tional refinements in order to adeguately represert the 
effect of switching transients on the transistor ani diode 
€lenents. 


BB. HALL EFFECT SENSCR FEEDBACK 


As was mentioned in the introduction, the switching 
action cf the commutation circuitry is based upcen rotor 
position feedback frem Hall effect sensors. In ercer to 
represent electronic switching and postional feedback, as 
well as the three phase winding configuration, it was neces- 
Sacvy te modify the tlock diagram of Figure 2.1 to that in 
Fiyure 4.3 

Comparing Figures 2.1 and 4.3, it is seen that the 
chief difference is in the switching logic block and the 
addition of two more transfer function blocks to acccunt for 
the additional windirgs. To better understand the logic 
relationship between the sensors and the current switches, 
the truth table for both counterclockwise and clockwise 
rotation is presented in Table II. 
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Figure 4.3 Equivalent Circuit With Electronic Commutation. 


A brief review of the action afa Hall effect senscr is 
in order at this point. Consider a conducting material with 
a constant current flow. When there is no external magnetic 
field present, the current flows undisturbed in a straight 
line. But when a magnetic field is applied, a deflecting 
force, the Lorentz force given by Il x B (where I ig the 
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TABLE II 
Sensor and Switching Logic 


RPS B RPS C PH A PH 5 PE C 
1 0 high low oyer 
0 8) high open iow 
Q 1 open high icw 
0 1 low high Oren 
1 1 low Open high 
1 0 open low high 

COUNTERCLOCKWISE 
0 1 low high oer 
0 1 oper high lcw 
0 0 high opén low 
1 0) high low oren 
1 ) open low high 
1 1 low open hich 
CLOCKWISE 


1 
magnetic field) causes the current to bend ané electrons to 


is the length of the conductor and8 is the 


current, 


pile up on cne side and positive charges, or holes, to dc 
the same on the other side. Thus a transverse Hall fpeten- 
tial difference, Vh, develops across the conductor (sez 
Figure 4.4). 

Figure 4.5 shows that i: the presence of a Sinuscical 
tagnetic fieli produced by a rotating magnet, such as cne of 
the pole-pairs of a multiple-pole motor, the induced Hall 
voltage takes cn an ideal square wave shape. 

The logic circuitry of the sensor is configured in such 
away as to output a logic level 1 for positive Hall voit- 
ages and a 0 for negative voltages. Briefly reviewing the 
acticn of the motor: as the shaft turns, the rotor position 
sensors send this information to the switching logic tlock 
which then decides which phases are to be macGe high, lcw or 
open in accordance with Table Il. Equivalently, the 
voltage, En, is then applied to the two appropriate windings 
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Figure 4.4 The Hall Effect 


and zero voltage is applied to the third winding in agree- 
ment with the assumption that current is steady in one leg 


while the commutation in the other two legs is eyuivalent to 


a continuous current in ne leg. As before, each voltage 
then produces a current. Because this is assumed to be a 
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Figure 4.5 Sensor Logic Based on a Hall Effect Senscr 


linear system, the principle of superposition is apflicatle 
and, thus, the torgues developed by current flow through 
separate windings car be considered independently and then 
summed to give a total motor torque. Beyond this pcint, the 
trushless motor comgenents are ilentical to its Erush-tyfe 
counterpart and behave similarly. 
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C. MCDEL REVISION 


CSHF modeling of this version of the brushless dc rotor 
called for two modifications to the brush-type aodel. 
Essentially, a CSME procedure function was constructed to 
act as a subroutine that embodies the logic of the switching 
action. The logic recognizes that if the motor shaft is 
within scme degree rarge, for example 30 to 60 degrees, then 
it sets the appropriate rotor position sensors and erergizes 
the switches for the appropriate two windings and deener- 
gizes the switch for the third. In effect, only the forcing 
voltages are passed to the main program, whereas switch an? 
sensor pesitions are available for output (see Appendix A). 
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A. AIR GAP FLUX 


The steady state behavior of both a brush-type and a 
-trushless dc motor has been studied under the fundamental 
assumpticn that the rack emf generated is a product of tae 
shaft velocity times a constant, Kb (see Fguation Z.z). 
Likewise, the motor torgue is equated to the product of the 
current times a constant, Kt (see Equation 2.12). But as 
Eguation 2.11 reveals in the case of the torque constant, 
Kt, the torque is in fact the product of some other constant 
cf£ properticnality, Kt1, times the air gap flux, $. The 
same is true for the tack emf constant, Kb; that is, Kb is 
the product of some proportionality constant, Xb1, times the 
air gap flux, 4. Ir both cases, the flux, ¢, has beer 
assumed tc be constant which is in agreement with the 
general frractice when modeling permanent magnet dc motors. 
In point of fact, the flux is not constant but rather is 
distributed sinusoidally, generally being the sum of a siru- 
soid and one or more of its harmonics. How the motor is 
physically constructed determines this flux distributicn. 
Two of the chief factors are the structure of the magnetic 
fole faces and, most importantly, how the windings are 
distributed in tne armature. The latter includes such 
factors as the number of slots, the number of coils fer slot 
and the spacing of the slots [Ref. 5]. The careful desijner 
distriltutes the windings so as to minimize the effects o£ 
harmonics. Regardless of their constituents, these 
composite sinusoids have an average value, of course, ani so 
in practice it is this averajye value that teceres the 
constant flux tera, ¢. 
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Be. FLUX AS AN AVERAGE VALUE 


A logical next step in the design of a amore realistic 
medel was the determination of an average value for the 
varible flux, davg, and the corresponding protortionality 
constant, KKi such tkat the ratio of davg and KX1 equals cane 
and the relationshi;ys of eguations 2.2 ani 2.5 are 
preserved. The equations for back emf, eb(t), and develorzei 


torgue, tm(t), now become 


eb(t) = KEp * wm(t) (egn 5.1) 


tm (t) Ktp * i(t) (egn 5.2) 


where the back emf ccnstant, kbp, and the torgue constant, 


Ktp, are now expressed as 


Kbp 


Kb * (avg/KK1) (egn 5.3) 


Kt * (avg/KK1) (egn 5.4) 


Xtp 


Assuming as a first approximation that the air gap flux 
varies as a simple sinusoid of unit magnitude, the first 
step toward finding the total system flux, avg, was to 
determine the relaticnship of the individual flux patterns 
in each of the three phase windings. Analysis of a four- 
role magnet rotating under three windings separated by 12) 


degrees (electrical) revealed the following relaticnshirs 


da = sin(26 + 1/6) (egn 5.5) 

db = sin(20 + 9176) (egn 5.6) 

dc = sin(2e + 5T/6) (egn 5.7) 
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where ¢a, @5. and ¢c are the flux contributions of the 
respective phases. Feturningjg to an eariz 

the circuit can be approximate] as two 

ignoring the third, then the average 

approximated froma the sum of tne flux develop 

ings at a time. This approximation is 

the fact that the flux in the decay le; ias, 

average value of zero during every conuutation 

From the switching seyguence of Figure Gt, 

conducting iegs were identified for eaca interval and, roa 
Equations 5.5 through 5.7, their respective flux values 
calculated and summed. For every interval, the resulting 
fiux waveshape was identical and had an averaje value of 
davj = 1.631 (see Figure 5.1). 

Since KK1 merely normalizes davg, its value was easily 
determined to be KK1 = 0.613. The appropriate substitutiors 
were then madie in the CSMP model (See Appendi A for 
Revision Two). Since this amounted to nothing more than a 
rearrangemert of corstants, the moiel'’s output behavior 
remained unchanged. At this stage, the nodel simulates tae 
behavior of a brushless de motor being electronically ccnau- 
tated and exhibiting constant, steady state developed tercue 
and motor speed due to a fixed, averaye value air gap flux. 


Ce. SINUSOIDAL FLUX 


The final stage of this thesis was to examine the 


steady state behavior of the trushless dc motor under vari- 
able air gap flux conditions. As a first approxitaticn, the 
flux was mojeled as a simpie sinusoid an? the flux relaticn- 
skips are those given in eyuations 5.5 tarough $.7 as 
before, only the £lux relatinj to tne two coniucting legs 
was considered. This composite variable lux, fvl, is shown 
in Figure 5.1 Eguations 5.3 and 5.4 nad to ke modified to 


account for the different value of flax. 
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Figure 5.1 Composite Plux Variation for Two Windings 


Kbpp Kb * ($v 1/KK2) (egn 5.3) 


Ktpr Kt * (dv1/KK2) (eqn 5.9) 


Recall tkat the objective was to keep the values of Kkpr and 


Ktpp aS near as possible to those of Kb and Kt, respec- 


tively. In order to roughly normalize ¢$v1, the constant 
value was determined as KK2 = 1.655. Since ?v1 was now a 
variakle guantity, it should be recognized that Kbprp and 


Ktpp were no longer strictly constants and were better 


expressed as the follcwing functions. 


Kbpp = Kb ($v1) (eqn 5. 10) 


39 


Ktpe = Kt (v1) (eqn 5.11) 


The eguaticns for back emf, eb(t), and developed torque, 


tm(t), now became 


— eb (t, dv1) 


Kbpp * wm (t) (egn 5. 12) 


tm(t,dv1) Ktpp * i(t) (egn S. 13) 
Fecause the total flux was no longer a constant, 
old average value, it was expected that variation in developed 
torgue (and, conseguently, in motor speed) would result. 
Furthermcre, as the system bleck diagram of Figure 2.1 indi- 
cates, torque is ccnverted by the motor transfer function 
ad into a retational speed. This is a process that is equiva- 
lent to inertia filtering of the toryue ripple that resuits 
from variable flux. For the purposes of this thesis, riprle 
is defined as the ratio of the amount of variation from the 
Cc Maximum to the maximum itself. Having made the aprrorriate 
substitutions in the model (see Appendix A for Revision 
: Three), the program was run under no-load conditiors. it. 
Ge was ckserved that there was, in fact, ripple in both the 
develcped torgue and motor speed and that some inertia 
filtering did occur. The next logical step was to icac¢ the 
motor and study the effect on torque and speed frifple. 
Table III indicates that under no-load conditions where 


speed was highest, the percentage of ripple in the outfut 
speed was at a minimum but that torgue ripple was at a 


maxinbun. 
It can be seen that as the load increased the amcunt of 


| 
ripple in the develored torgue decreased from about 84% to “wa 
1 
8%. A Closer examination of the torque behavior under : 


varying loads revealed that, even as the aotor torgue neces- 


sarily grew with increasing loads, in every case there was 
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TABLE III 
Torgue and Speed Ripple Due to Sinusoidal Flux 


MOTCR MOTOR % RIPPLE 

IOAD (oz-in) SPEED (EFM) TORGUE SPEED 

0.0 2557 33.9 1.1 

32.0 2087 21.1 1.6 

€4.6 1617 12.2 2.6 

96.0 1147 7.8 4.6 
approximately a constant 8-9 oz-in variation. This varia- 


tion can be attributed primarily to the flux variaticn ir 
Ktpp and to a lesser degree to the variation in i(t) caused 
ky the variation in tack emf. Thus tne magnitude of torjue 
variation remained akcut the same for all loads though the 
percentage ripple arpears to indicate otherwise. In the 
case cf motor speed, the amount of ripple appeared to 
increase slightly with increasing load and thus decreasing 
speed. In this instance, the magnitude of variation did 
increase a small amount with decreasinj speed - from 2° rpn 
at no-load to 52 rpm at 96 oz~ins - and suggests that speed 
tipple, at least, is somewhat depenient on load. This is 
due in part to the fact that as the motor siows, there i 
less rotor momentum and thus less inertia filtering. Ta 
addition, at slower sreeds there are less frictional effects 
which at higher speeds also tend to filter out ripple. in 
general, it makes intuitive sense that at higher sreeds 
freguency variations are harder to discern than at slower 
speeds. 


D. HARACNIC FLUX 


Having observed the effects of simple sinusoidal flux 
variation, the next step was to replace the single sinusoid 
with a flux composed of the sum of a sinusoid and its 


harmonics which is mcre representative of the air gap flux 


patterns in an actual motor. The back emi of a typical 


commercial Erushless dc motor is given in Figure 5.2 where 
the voltage between two terminals was measured with the 
motor rotating at 1200 rpm. 

A harmonics analysis of this waveshape shows that the 
principle harmonic was the fifth and that the waveshape was 
closely described by the expression 


3.0 sin(6) + 0.59 sin(5é) (egn 5. 14) 
When sukstituted into eguations 5.5 thru 5.7, and the 


different angular speed taken into account, the fcollewirg 
relationships for flux result 


da = 3.0*sin (20+176) + 0.59*sin (100+5 17/6) (egn 5. 15) 
bb = 3.0%sin(29+9 1/6) + 0.59#sin(109+9 1/6) (eqn 5.16) a 
Oc = 3.0%sin(20+5 T/E) + 0.59*sin(108+ M6) (eqn 5.17) i 


The same procedtres used in the simple sinusoid nodel 
were then applied to the harmonics case where the compesite 
flux changes to $v2, a flux factor that varied to a slightiy 
greater dJegree than did dvi. Equations 5.8 and 5.9 hai to 
ke mcdified as follows 


Kbppp = Kb *® (Ov2/XR3) (egn 5.13) 
Ktpep = Kt * (QV¥2/KK3) (eqn 5.19) 
where KK3 approximately normalizes the flux and had the 


value KK3 = 5.38. As before, Kbppp and Ktppp varied with 
the flux and so the new equations for back emf and develoreda 


torgue becane 
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Figure 5.2 Back EMF Waveform Due to Harmonic Flux 
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eb (t, v2) Kb (Ov2) * wm(t) (egn 5.20) 


eb(t,$v2) Kbppp * wm(t) (eqn 5.21) 
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Kt (Ov2) * i(t) (egn 5.22) 


tm (t, v2) 


tm (t,dv2) Ktppp * i(t) (egn 5.23) 

The appropriate substitutions were once more made in the 
model (see Appendix A for Revision Four) andthe modjei was 
run under various loads. As anticipated, the variable flux, 
¢v2, produced ripple in the developed torgue and motor 
speed. Table IV shews the resultant values for various 


loads. 
TABLE IV 

Torgue and Speed Ripple Due to Harmonic Flux 
MOTCR MOTOR & RIPPLE 
LOAD (oz-in) SPEED (rpm) TORQUE SPEED 
0.0 2557 94.6 2.4 
32.0 2087 38.5 3.7 
64.0 1617 22.5 5.6 
96.0 1147 14.7 9.8 


Comparing Tables III and IV, it can be seen that the values 
for torgue and speed ripple are of Similar magnitudes. The 
somewhat higher values for torgue are due to the fact that 
in the $v2 case there was a 16-21 0z-in variation as offosed 
to an average 8-9 cz-in variation in the Ovl case. The 
speed ripple behavior was much like the ¢v1 case with 
proportionate increases resulting from the increased torgue 
ripple. Thus, for both torgue and speed ripple due to v2, 
the model performed as expected; that is, slightly higher 
values of ripple occurred which were clearly attributable to 
the increased variability of the flux pattern of ccmpcsite 
Sinusoidal harmonics. 

An important observation made during runs of the 6v1 
model came as a consequence of the motor speed variability, 
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particularly when tke motor was unloadied. Since no-load 
speed was a function of Fb as given in equation 3.4, the 
vaciakle Krtpp had the effect of increasing the upper-end 
speed beyond the 2557 rpm baseline given in Figure 3.1 This 
increased magnitude ccupled with the swing of speed values 
created an undesirable effect in the unloaded conditior. In 
effect, at certain points in time the additional speed 
produced tack emf values that exceeded the input voltage 
which in turn resulted in negative current and developed 
torque values. Though these were transient negative values 
and had small effect on the steady state operation of the 
motor, still their presence was obviously unrealistic. The 
soluticn to returning the motor to all positive values of 
current and torque was to increase the fricticn filtering 
effect by adjustment of the viscous friction coefficiert, 
Bm. It will be recalled from Chapter Three that a value for 
the Bm of the commercial motor being modelled was not avail- 
able and that for the basic model a value for 8m was derived 
from the curve fitting process. Thus, it seezed reasorable 
and pérmissable to refeat the procedure here. 
Since it makes sense that increasing the friction within 
a system will slow the system down, in this case increasing 
the value of Bm did just that. By increasing Bm fren 
0.00015 to 0.045 oz-inyrad/s, the upper-end of the motor 
speed was reduced to a value that eliminated the undesirable 
negative values. This caused the no-load speed tc te 
lowered slightly belcw the desired value of 2557 rpns. 
Recalling that no-lcad speec€ is fixed by the back enf 
constant, Kb, this parameter was adjusted from 0.1120 to 
0.1089 velt/rad/s to kring the back up to speed. Since this 
adjustment was less than 3% and the manufacturer specifica- 
tions allowed a 10% +/- measurement error margin, this 
adjustment to Kb was pernissable. Of course, a propor- 
tionate adjustment to the torgue constant, Kt, was also made 
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since Kt and Kb must maintain a constart relationship as 

explained in Chapter Three. The effect cn gotor speed _ 

Tripple was neglible. Reviewing the $v2 model at this roint 

revealed that similar behavior in the no-load state occurred 

but cn a larger scale due to the greater variability of v2. : 

A Similar set of adjustments was necessary to elininate the os 
negative current and torque values in the unloaded cendi- 
tion. It reguired an adjustment of Bm from 0.00015 to 0.045 
oz-inyrad/s in order to ensure that all values were jfosi- 
tive. In addition, the back erf constant, Kb, was adijusted 
to 0.1192 voltyrad/s and a proportionate acjustment to Kt 
was made. Again, all adjustments were within tolerances. 

Of course, it was earlier established in Eguatior 3.7 
that changes to Kb and Kt would change the slope of tke 
speed-torgue curve. Thus, in both cases the slope was 
returned to its original value by adjusting the motor resis- 
tance, R. In the case of sinusoidal variable flux, & was 
reduced from 2.74 to z.70 ohms. For harmoric variatle flux, 
R was increased from 2.74 to 3.0 ohms. Again, ktoth adjust- 
ments were within the 10% +¢/- allowable deviation. It was 
also cbserved that since both Kb and Kt no ionger remained 
constant but rather varied within small ranges, the slope 
was not constant and varied accordingly in both the $v1 ana 
¢v2 models with the greatest change occurring wher Kk ard xt 
Simultaneously reached the minimum values of their respec- 
tive ranges. Bearing in mind that these deviations in slope 
are fairly transient, their average values are most impcr- 
tant and are the apfrfroximate values from which the final 
adjustments to R were based. 


46 


Tae . . . . “. ee a Se - ae = oy oe ie = oy - le . os @ « ae a > ra Ps ial a . . aa 

Loe ne eee erate Toe Bose oe ve aieiare Seige oa et a he ee ee Se Seg a gi Sig ia Stele cate hartge aa nin re om 
oe age eee re eto ee 5 ag fae a ee Se ea Seed we . 
an ne a ear nr Ye en a atte eta as wee ot ea a a a 


TO 4 e = ~ wo oe ee Je See: ek 


a VI. SUMMARY 


A. REMARKS AND CONCLUSIONS 


CSMP language is a convenient tool for modéling a 
Erushless DC mQqtor. Positional seasor feedback and 
switching logic being functions of time, the action cf €iec- 
tronic commutation can be nicely Simulated in CS*E&. Ihe 
model can be fairly easily modified in order to study cther 


& | 


motor configurations. For example, the three-phase star 
configuration used in this thesis could easily be changed to 
a grounded neutral, point Din figure 4.1, and the same 
types of analysis performed. In future stujies where the 
complexity of the mcdel is increased in oréer to include 
more design detail, the advantage of CSMP will become even 
more apparent. 

In the course of simulating the effects of different 
air gap flux variaticns, it appeared necessary to make @irnor 
parameter adjustments in order to offset the feedoacs 
effects cf torque and speed ripple. Within the scope of 
this thesis, adjusting the steady state performarce «as 
reasonable and justifiable. In a larger context, however, 
various control systems, particularly torgue generation 
contrcllers, are specifically desiyned to reduce tccyue 
tipple and its effects. So, in at least one sense, thes2 
adjustments were somewhat artificial. To the extent that 
parameter adjustment is necessary and desirable, it must be 
understood that this is not a simple procedure since adjust- 
ment of one parameter aimost always requires adjustment of 


one cr more others. These adjustments in turn necessitate 


readjustment of the criginal parameter. The process is thus 


an iterative one. This model permits this kind cf parareter 
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adjustment, but it is a long and tedious process that would 


ke ketter performed Ey some scrt of optimization algoritau. 


BE. RECCMMENDATIONS FCR FPOTURE STUDY 


The simulation cf ripple effects is particulary rele- 
ij vant to the control engineer who must design the control 
circuitry and power ccnditioner to meet specific performance 
requirements. This points to several areas fer future 
study. One area is the contrel of torque jeneration and the 
several principle configurations used in brushless dC 
motors. Two commonly used methods are the sinusoidal torgue 
generation principle and the trapezoidal torque scheme. 
(Ref. 5]. In addition to torgue generation control, a 
second area for investigation is power control of brushless 
PC motors. In this thesis, output power was unregulated ani 
was studied simply as a function of a constant inyut 
voltage. In practice, however, power is produced ani 
controlled by varying the supply voltage. Tarious schezes 


Bi 


exist fer power control; amoung these are iinear transistoc 
contrel and pulse-width or pulse-frequency contrel. 
[Ref. 5}. 
= A very important area that this thesis did not address 
was that of the switching transients that occur as a conse- 
guence of commutation. During each Switching interval, each 
winding is an inductcr that stores energy which causes a 


Significant amount of current to flow at the instant of 
switching. This can have serious effects, especially if 
breakdown conditions exist amoung the commutation transis- me dg 
tors and other contrel circuit elements. In this same area, - 
Switching transients also occur within the controi circuitry 
itself. Thus, a thorcugh investigation and understanding o£ 
the behavior of power transistors, diodes and other control 
elements used with electronic commutation needs to be dene. 
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LISTING OF MODEL PROGRAMS 


A. BASIC PROTOTYPE MODEL 


The program in this section is a basic program that 
Simulates the action cf a standard brush-type Dc moter. mt 
is the prototype from which the brushless DC zoter modei ani 


its several revisions are derived. 
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B. REVISION ONE 


This is the first revision of the basic prototyze. This 
model simulates the switching logic of a brushless [C rotor 
wnich is based on feedback from Hall effect sensors. an 


addition, the windings are treated indepenientl 


ne 


a 
contributions to developed torgue are superposeé. Since the 
superfesition results in twice as much current flew asa 
single iumped coil, the total current is halved to retain 


the same overall moter behavior as in the prototype. 
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C. REVISION TWO 


In the following program, the motor's air gar flux is ae 
treated as the average value of the sum of the flux acting t 
in two windings as explained in Chapter Five. For cerven- aie 
ience of study, both windings are again lumped into a single a 
winding as in the frototype. This facilitates a closer a 
focussing on the effects of variable flux on motor Lehavicr. 
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De. REVISION THREE 


The purpose of this version oz the program is treat 
air cap flux as varying in simple sinusoidal fashion. 
total system flux is treated as the alyetraic sum of 
flux developed in the two active windings. 
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ECLRE PROVICES A SIMPLE MECHANISM FOR REVERSING THE MCTOR*S 
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E. REVISION FOUR 


This version of the model treats tne flux as varying 
according to the sum cf a sinusoid of fundamental frequency 
and its fifth harmonic as explained in Chapter Five. the 
total flux is again approximated as tne algebraic sur cf the 
flux developed in twe windings at a time. 
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SAMPLE OUTPUT 


The following is a sample of the type of information 
the mcdel can provide. In this particular case, the fourth 
revision of the basic prototype, the harmonic flux case, was 
used. For completeness, the motor was run at loads fro 0.9 
up through 96.0 oz-in, which was the range of Ilcads with 
which all analysis was performed. In addition to such motor 
variakles as speed, current and power, a sampling of the 
feedback from Hall effect sensors as well as the switching 


action of the power transistors is included. 
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